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Deformation of a flexible fiber in a viscous flow past an obstacle 
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We study the deformation and transport of elastic fibers in a viscous Hele-Shaw flow with curved 
streamlines. The variations of the global velocity and orientation of the fiber follow closely those of 
the local flow velocity. The ratios of the curvatures of the fibers by the corresponding curvatures of 
the streamlines reflect a balance between elastic and viscous forces: this ratio is shown experimen¬ 
tally to be determined by a dimensionless Sperm number Sp combining the characteristic parameters 
of the flow (transverse velocity gradient, viscosity, fiber diameter/cell gap ratio) and those of the 
fiber (diameter, effective length, Young’s modulus). For short fibers, the effective length is that of 
the fiber; for long ones, it is equal to the transverse characteristic length of the flow. For S p < 250, 
the ratio of the curvatures increases linearly with Sp ; For S p > 250, the fiber reaches the same 
curvature as the streamlines. 


I. INTRODUCTION 

The transport of flexible biological or man made fibers by a flow and their deformation is of interest in view 
of their potential applications in many different industrial fields. A first example is the paper industry in 
which the orientation and spatial distributions of the fibers in the flowing pulp must be controlled^. Fibers 
are also widely used by petroleum engineers^ to enhance the proppant transport capabilities of fracturing 
fluids, or to prevent the backflow of proppant. Other domains of application are civil engineering (special 
cements, structural reinforcement), textile engineering, bio engineering and medicine. 

In flows at low Reynolds numbers, the deformation of fibers reflects a balance between viscous and elastic 
stresses and plays often an important part: for instance, these effects are critical for understanding the 
dynamics of flexible biological filaments^—, the filtration and cross stream migration of macromolecules 
in small channels^— and the positioning of bio-fibers in porous medial. In all these applications, the 
prediction of the motion of the fibers is a difficult fundamental problem, particularly in complex geometries 
and inhomogeneous flows where constrictions and obstacles are present: these create complex interactions 
between the flexible objects, the flowing fluids and the solid walls. 

Our interest in these problems is focused on the injectability and transport of long flexible fibers through 
fractures in reservoir rocks. In a previous workH, we studied fiber transport by a flow in single narrow model 
rough fractures and of its dependence on the fluid velocity, the flexibility of the fiber and the configuration 
of the aperture field. 

The present work is devoted to a quantitative study of the influence of the parameters in the configuration 
where the fiber passes in the vicinity of contact areas or close to regions of low aperture. This situation is 
modeled experimentally by means of a Hele-Shaw cell of constant aperture and of width varying with the 
distance: this allows one to determine the variations along the flow of the orientation and deformation of the 
flow lines and of a fiber placed in the flow. More specifically, we seek to relate the variations of the curvature 
of the fiber to those of the streamlines, to the velocity gradients and to the geometrical and mechanical 
properties of the fiber. 

Most past studies considered the simpler case of flexible fibers in a non confined parallel shear flowi^r— 
of shear rate ( dv / dn) (v is the modulus of the local flow velocity and the derivative is taken is the direction 
perpendicular to the flow lines). In such a configuration, the fiber may rotate about its axis and buckles 
above a threshold value of the Sperm number which characterizes the relative magnitudes of the viscous and 
elastic forces, and is equal to: 

p,(dv/dn)£ 4 ' 
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fj, is the viscosity, ( dv/dn ) is the shear rate along the fiber, E its Young’s modulus, / its area moment of 
inertia (/ = 7 tD 4 /64 for a circular cross section of diameter D) and £ a characteristic dimension (for instance 
the length Lf of the fiber). Other studies dealt with the cross stream migration of a flexible fiber in a parallel 
Poiseuille flow which is also strongly influenced by buckling^—. Finally, buckling is also observed above a 
critical value Sp ~ I20i£) of Sp for fibers approaching a stagnation pointier—. Numerically, slender-body 
theory has been used to model the behavior of an elastic fiber (with one clamped end) in the vicinity of a 
corner in a curved channel^. In this latter work, the variation with time and the final shape of the filament 
were analyzed in different corner geometries as a function of Sp. Recently, Wexler et al M 2 considered a fiber 
with one clamped end and initially transverse to a flow confined between two plates (like in the present 
case). In all cases, bending results from the combined effects of the hydrodynamic and viscous forces on the 
fibers: yet, in Refs. [2l| and [ 22 ], their resultant is balanced by a force on the attachment point while, in the 
present case of a free fiber, this resultant must be zero, leading to very different shapes of the fibers. 

In the following, we characterize first the flow field in the Hele-Shaw cell of varying width by means of dye 
tracer measurements and of 3 D numerical simulations: we study in particular the variations of the velocity 
along the streamlines. Then, the local velocity and mean orientation of the fibers are compared to those of 
the fluid velocity at the same location. The curvature Kf of the fibers is characterized by its ratio k/ /k s 
to that of the streamlines: the variations of the extremal values of this ratio on different streamlines are 
studied as a function of the diameter D and length Lf of the fiber, of its initial transverse location and of 
the flow velocity. The definition of the Sperm number Sp of the problem allowing one to take into account 
the influence of these different parameters is then discussed and the experimental dependence of Kf/n s on 
S p is analyzed. 


II. EXPERIMENTAL PROCEDURE 


The experimental setup consists of a Hele Shaw cell placed vertically. Its height and gap are respectively 
L = 280 and H = 2.2 ± 0.1mm (Fig. [[]). The upper part of the cell has a Y-shaped section; the top is 
attached to a fluid bath with a slit at the bottom allowing for the flow of the fluid and the injection of the 
fibers. The fiber is injected in a funnel-shaped device which permits to select conveniently the location at 
which the fiber is injected. A rounded off step is then cut out of a transparent mylar sheet of thickness equal 
to H and inserted in the Hele Shaw cell. The mylar is tightly fitted so that no flow is observable between its 
surface and the walls of the cell. The height AW of the step parallel to the y axis is 30 mm and it is located 
in the middle of the cell; the width W (x) of the flow section in the cell is then 70 mm in the inlet section 
and is reduced to 40 mm at the outlet as shown in the figure. The axis y = 0 corresponds to the boundary 
of the mylar sheet at the inlet. 



FIG. 1. Experimental setup. 


Using a continuous rather than an abrupt variation of the height of the step boundary with the distance 
first minimizes unwanted solid-solid interactions between the fiber and this boundary. Moreover, a continuous 
variation of the width of the flow channel with distance is closer to the variation encountered in a real 
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fractured rock. Finally, putting uniform flow sections before and after the step allows one to detect possible 
irreversible rotations and deformations of the fiber. 

A gear pump produces a downward flow; two Newtonian fluids are used: pure (MiliQ) water, of density 
p = 997 ± 1 kg.m~ 3 and a viscosity p ~ 1 mPa.s, and a solution of 70% glycerol by weight in water, 
of density p = 1178 ± 1 kg.m~ 3 and viscosity p « 15mPa.s (27°C). The flow rate Q is in the range 
10 4 < Q < 2.3 x 10 4 mm 3 .s _1 , which corresponds to a characteristic velocity 95 < U < 230 mm.s -1 ; in this 
work, U is the maximum of the Poiseuille velocity profile (half way between the walls) in the inlet parallel flow 
section upstream of the step. Downstream of the step, the velocity is higher by a factor 1.75. For the water- 
glycerol solution, the corresponding Reynolds number Re = UH/v is in the range 16 < Re = UH/v < 40. 
For water, one has respectively 210 < Re < 500. 

The flexible silicone fibers are made using a technique developed in the laboratory. Silicone is forced out 
of a syringe at an adjustable flow rate by means of a syringe pump. A second pump with a plate attached 
to the top is placed below the end of the first one so that the plate moves perpendicular to the flow. The 
injection leaves on the top of the plate thin flexible silicone fibers of constant diameter. The respective values 
of the speeds of the silicone at the outlet of the syringe pump and of the moving plate are very important. 
In order to obtain fibers of constant diameter, the velocity of the outflow must match approximately that of 
the plate. If the speeds are not well matched, either the stream of silicone buckles or it breaks^. Fibers of 
different diameters can be made by using nozzles of different sizes: in the present work, nozzles of diameter 
300 — 600 pm were used. The average D of the diameter over its orientation in the section and the fiber 
length is constant within ±0.015 mm from a fiber to another; the diameter is constant along the length for a 
given orientation within ±0.05 mm (D = 0.49mm) and ±0.025 mm (D = 0.49mm). In a given section, the 
variation of the diameter with its orientation is of the order of ±0.15 mm for D = 0.91 mm and ±0.02 mm for 
D = 0.49 mm. Large fibers become indeed elliptical because they sag under gravity before they polymerize: 
this effect is prevented by surface tension for D < 0.7 mm. 

The lengths of the fibers are Lf = 15, 30 and 60 mm and, although they are lighter than water (pf = 
796 kg.m -3 ), the tests showed that gravity has little effect on the experimental results. In order to evaluate 
the stability of their properties, the silicone fibers were left overnight in different fluids (water, and water 
glycerol solutions): no swelling occurred. 

In the interpretation of the present experiments, the bending stiffness El of the fiber plays a major role. 
In order to measure E /, the fiber is allowed to bend under its own weight with one of its ends clamped in a 
chuck-like device which keeps it locally horizontal; the fiber is placed in front of a dark background in order 
to achieve a high optical contrast. The process is repeated for several fiber lengths: 5 < Lf <80 mm. The 
values of Sx, i.e. the vertical deflection of the tip of the fiber with respect to its position or a zero external 
load, and of Lf are determined with an uncertainty of ±20 pm by means of an image analysis technique 
developed by Semin et al^ 4 . In the present case of a rod bent under its own weight and clamped at one end, 
the theoretical expression of Sx is mi gL“f/{8EI) where meg is the force per unit length (i.e. the weight 
per unit length). The data were analyzed by plotting the deflection Sx as a function of the fourth power L 4 
of the length. The relationship between these two quantities is linear for Sx/L < 0.1; this is the limit of the 
small deflection approximation. The slope of this plot gives the value of El (See Tab]]]). 


Fiber 

Material 

D (mm) 

P (kg.m 3 ) 

El x 10 B (N.nF) 

FI 

Silicone 

0.49 ±0.015 

0.8 

0.28 ±0.04 

F2 

Silicone 

0.91 ±0.015 

0.8 

0.37 ±0.04 


TABLE I. Physical and mechanical characteristics of the fibers ( D , p and El are respectively the diameter (averaged 
over the length of the fibers and the orientation in a section), the density and the bending stiffness of the fiber). 

Experimentally, after the flow is established, the fibers are inserted into the upper inlet with one of the 
ends held by pliers. Three insertion points corresponding to different distances ymiet from the boundary of 
the step at the inlet have been selected (see Fig. (2^). Once the fiber is aligned with the flow, it is released 
and digital images of its subsequent motion in the cell are recorded. In order to improve the contrast, the 
model is illuminated from behind by a plane light panel, and images are obtained by means of a digital 
camera at a rate of 30fps and with a resolution of 1024 x 768 pixels (with lpixel = 0.22 mm). Each picture 
is then processed digitally in order to determine the instantaneous profile of the fibers. 

The velocity of the fiber is determined from the displacement of its geometrical center between two 
consecutive pictures. The deformation of the fiber is characterized quantitatively by adjusting locally the 
profile by a polynomial function in order to determine the local slope and the curvature. Only experiments 
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in which fibers are relatively straight and parallel to the flow before reaching the step were considered in 
the analysis. 

In parallel with these experiments, three dimensional numerical simulations solving the full Navier-Stokes 
equation by means of the FLUENT™ package have also been carried out. These simulations provide both the 
components of the velocity and their derivatives with respect to the coordinates x, y and z. The curvature 
of the streamlines and the derivative dv/dn of the velocity with respect to distance perpendicular to these 
lines is then computed by combining the values of these derivatives. 


III. EXPERIMENTAL RESULTS 
A. Characterization of the flow field 



FIG. 2. (a) Visualization of three streamlines by the injection of a dye (continuous lines). Gray region: step obstacle; 
dotted lines: streamlines computed numerically for the same injection points, (b) Numerical (lines) and experimental 
(symbols) variation of the fluid velocity v along streamlines corresponding to three injection points located at x = 0 
and yiniet = 4.5 ± 1.5 mm ((o), dashed line), 13.5 ± 2.0mm (A, continuous line) and 21.5 ± 2.0mm (□, dotted line). 
The characteristic velocity U upstream of the obstacle is equal to: 230 mm.s -1 . Note that the horizontal and vertical 
length scales in Fig. a are different. 


Prior to the injection of the fibers, we have studied the flow velocity variations along the streamlines 
starting at each of the three selected injection points. Experimentally, a stationary flow is first established in 
the cell and a stream of dye is then injected at each location (Fig.[2]i). The progression of the dye is recorded 
and analyzed using the procedure described in the previous section. The tip of the dye streak moves at the 
maximum of the local Poiseuille profile in the gap: this allows one to measure quantitatively the variation 
of this maximum along the corresponding streamline. The results obtained from these measurements are 
then compared to those of the 3 D numerical simulations. 

Figs. [2]a-b show that the experimental and numerical results are in good agreement for both the geometry 
of the streamlines (a) and the variations of the fluid velocity along them (b). For the most distant streamlines 
(dotted line and squares in Fig. [5]), the velocity increases continuously with distance along x in direct relation 
with the decrease of the flow section. The results are different for the two other streamlines: the velocity 
first drops as the fluid reaches the vicinity of the step, then it increases sharply, reaches a maximum right 
after the step and decreases finally towards a constant value. Localized buckling may therefore occur in 
these small regions of decreasing velocity. 

Fig. [5] corresponds to an upstream velocity U = 230mm.s _1 (in the mid plane between the walls): 
these numerical simulations were also performed for a velocity 10 times lower. As expected from the Hele 
Shaw approximation, the streamlines of interest at both velocities (and their curvature) are identical within 
experimental error while the velocity components are everywhere proportional to U. 

Actually, the deformations of the fibers are not due to the v velocity itself but to its spatial variations 
(see Section HUDl) . The variations with the distance x of the normalized gradient components (dv/ds)/U 
and (dv/dn)/U parallel and transverse to the velocity are plotted in Figs. [3Ji-b for the same streamlines 
as in Fig. [5}j. The parallel component (dv/ds)/U has its highest positive value (inducing an extensional 
force) in the region of maximum slope of the streamlines (Fig. [3Jr); it is instead negative (and produces a 
compressional force which may induce buckling) at the very beginning and end of the steps. The transverse 
component ( dv/dn)/U which will be seen to induce deformations of the fiber (see Sec. IIII Cl) goes instead 
to zero when the slope of the streamlines is highest; ( dv/dn)/U reaches its maximum positive and negative 
values in the regions of maximum curvature (Fig. [3J 3 ) - The width of these extrema is of the order of 30 mm. 
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FIG. 3. Variations as a function of the distance x of the normalized parallel and transverse components (dv/ds)/U 
(a) and (dv/dn)/U (b) of the gradient of flow velocity v computed numerically for the same streamlines as in Fig. [2j> 
(the different types of lines have the same meaning as in this latter figure). 


Quantitatively, Tab. HT1 lists the minimal (negative) and maximal values (as determined from the numerical 
simulations) of the curvature of the line together with the transverse normalized velocity gradient ( dv/dn)/U 
at the same locations; as mentioned above, these values are independent of U. At the distances x = 90 and 
130 mm corresponding to the minimal and maximal velocities along the streamline closest to the obstacle, one 
finds respectively: ( dv/dn ) ~ 0.01517 and 0.0417. The experimental values are very close to the numerical 
ones. 


Init. location 

Max. curvature (upstr.) 

Min. curvature (downstr.) 

yiniet 

max 

(1 /U)(dv/dn) 

mm 

ni s 

(l/U)(dv/dn) 

4.5 

0.0144 

-0.014 

-0.0212 

0.039 

13.5 

0.0109 

-0.006 

-0.011 

0.0222 

21.5 

0.0077 

-0.0034 

-0.0066 

0.011 


TABLE II. Transverse coordinate yi n iet (in mm) at inlet, maximum (k™ ox ) and minimum (ft!) 1 *") curvatures and 
velocity gradient (1/U)(dv/dn) in mm -1 at the corresponding locations for the three streamlines shown in Fig. [5] 


B. Mean velocity and orientation of the fiber 


50 

Bf,6 s 

(°) 

25 


0 

0 100 x (mm) 

FIG. 4. (a) Experimental variations as a function of the distance x of the flow velocity v measured by the injection 
of dye (lines) and of the velocity Vf of the center of mass of a 30 mm long fiber of diameter D = 0.49 mm (symbols) 
along streamlines corresponding to the three same injection points as in Fig. [2] (b) Variation as a function of the 
distance x of the local angles with the axis x of the flow velocity (lines) and of the fiber (symbols). Fluid velocity 
upstream of the obstacles: U = 230 mm.s - . The meaning of the symbols and of the different types of lines is the 
same as in Fig. 

We have studied the velocity Vf of the center of mass of fibers of diameters D = 0.49 and 0.91 mm 
corresponding to confinement ratios D/H = 0.22 and 0.41 and initially straight and parallel to the flow. 
In the parallel flow region upstream of the step where v = U, one finds that: Vf = (0.97 ± 0.05)17 and 
Vf = (0.87 ± 0.1)17. The ratio Vf/U only depends on the diameter of the fibers and not of their length. 
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Moreover, from Fig.[5Ji, the velocity Vf remains close to v at all distances along the streamline and not only 
in its parallel part. These results agree with an experimental and numerical study by other authors^ of the 
transport of rigid fibers by a Poiseuille flow between parallel flat walls. For the same ratios D/H = 0.22 and 
0.41, these authors predict that Vf is respectively only 4% and 10% lower than v. Vf is, in addition, the 
same for fibers parallel and perpendicular to the flow. This small difference between Vf and v, particularly 
for D/H = 0.22, implies that the friction forces due to the relative motion of the fiber and the walls do 
not play a major part in the phenomenon. The velocity of the fluid in the center part of the cell gap is, 
therefore, used as the relevant velocity of the problem. 

As the fiber moves along the streamline, it becomes tilted with respect to the vertical axis x as it moves 
over the step: the variations of the tilt angle 0 / with the distance x for the three streamlines are plotted in 
Fig. [4l together with those of the local tilt angle 9 S of the streamline at the same distance. The variations 
of Of and 9 S are very similar, particularly for the two streamlines farthest from the contour of the step. For 
that closest to the contour, instead, the maximum tilt angle of the fiber is slightly lower than that of the 
streamline; moreover, Of does not revert exactly to zero after the step. This suggests that a part of the fiber 
moves into the layer close to the contour where the Hele Shaw approximation of zero 2D vorticity in the 
(x, y) plane is not valid any more, resulting in a residual rotation. When the fiber happens to be injected 
particularly close to the obstacle, this rotation may be so large that the fiber touches the obstacle: such 
experiments are not included in the present interpretations. 


C. Deformations of fibers of constant length 



FIG. 5. Fibers (Lf = 30 mm, D = 0.49 mm) observed at two particular locations upstream (left) and downstream 
(right) of the step for three upstream velocities: U = 95mm.s _1 (black), 160 mm.s -1 (light gray) and 230mm.s -1 
(gray). Streamlines correspond to yMet = 4.5mm (dashed line) and yi n i e t = 21.5mm (dotted line). Fluid viscosity: 
fi = 15 mPa.s. 

In this section, we compare the deformations of fibers of same length Lf = 30 mm at different distances x 
on a given streamline and for three different values of the mean flow velocity. Fig. [5] displays superimposed 
images of fibers of same diameter D = 0.49 mm obtained at these different velocities in the regions of 
maximum positive (a) and negative (b) curvature of two different streamlines (these regions are located 
respectively on the upstream and downstream sides of the step). In case (b), increasing the fluid velocity 
strongly increases the curvature of the fiber which becomes of the order of that of the streamline: this is 
particularly visible on the streamline closest to the contour of the step. For the other, the curvature is 
always similar to that (weaker) of the other streamline. In case (a), the effect of the velocity is still clear 
close to the contour but much weaker on the other streamline. 

These informations are complemented by Fig. [6] displaying the variation with distance of the curvature 
Kf of the same fiber as it moves along the streamline corresponding to yiniet = 13.5 mm. This variation is 
compared for different velocities U to that of the local curvature k s of the streamline: Kf reaches extremal 
values at distances x close to those corresponding to the extrema and The minimum of Kf 

downstream of the obstacle is of the order of the minimum value of k s at all velocities while it is significantly 
lower and increases with U for the maximum upstream. Up to the end of the step, the curvatures of the fiber 
and of the streamline follow a similar trend of variation with the distance x, even though the amplitudes 
are generally different. Farther downstream and at the highest velocity, the curvature of the fiber displays 
a second maximum while that of the streamline returns to zero. This may be due to buckling induced by 
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FIG. 6. Curvature K/(mm _1 ) of a fiber (D = 0.49mm, Lf = 30mm) as a function of the coordinate x of its 
geometrical center for an initial location yiniet = 13.5 mm and for flow velocities: (O) U = 95 , (□) 168 and (o) 
230mm.s~ 1 . Continuous line: curvature of the corresponding streamline at the same distance. Fluid viscosity: 
fi = 15mPa.s. K ™ ax , K™ ln •. maximum and minimum curvatures of the streamlines. 


the locally negative value of the derivative dv/ds shown in Fig. [3Ji and resulting in a compressional stress 
on the fiber just downstream of the step. Figs. [5] and [G] both show therefore that the absolute curvature of 
the fiber increases with the local curvature of the streamlines but depends also of the local flow. 

We study now quantitatively the values of these extremal curvatures at different flow rates. Since the flow 
is laminar, the curvature of the fibers corresponds to a balance between elastic and viscous forces: the latter 
are induced by local velocity differences between different parts of the fiber and the fluid. These differences 
are small compared to the characteristic flow velocity U since the global fiber velocity is close to that of the 
fluid: the corresponding Reynolds number is also small and the viscous stresses are proportional to these 
velocity contrasts. The relation of the latter to the spatial variations of the velocity and the deformation 
of the fiber and the streamlines is discussed in Sec. IIII Dl but, as a first step, fi dv/dn will be chosen as the 
reference viscous stress: dv/dn is indeed strongly correlated to the variations of the curvature, as shown by 
Figs. [3] and [Gl 



FIG. 7. Dimensionless curvature, k f /k s , of fibers of length Lf = 30 mm as a function of the viscous stress fi(dv/dn). 
Diameter D = 0.49mm; (□): viscosity fj, = 15mPa.s, (^): g = 1 mPa.s. (o): D = 0.91mm, p = 15mPa.s . 

Fig. [7] displays variations of Kf/n s as a function of the viscous stress \xdv/dn ( Kf/n s is used as the 
parameter of interest since, from Fig. [6l one has always Kf < k s ). Significant curvatures kj/k s are only 
observed for the water-glycerol solution (viscous stresses are too low for water (* symbols)). As expected, 
for a given fidv/dn , the curvatures are also highest for the most flexible fiber (D = 0.49mm). In this latter 
case (□ symbols), k//k s increases at first linearly with p, dv/dn and tends then toward a limit of the order 






























of 1 (the fiber coincides locally with the streamline). This limit is not reached for the less flexible fibers (o 
symbols). 


D. Characteristic dimensionless number 

In order to collapse the data corresponding to different fibers onto a single master curve, it is necessary 
to replace as the horizontal scale fi(dv/dn ) by a combination including, in addition, the diameter and the 
length of the cylinder and the stiffness of its material. As mentioned above, the degree of bending of the 
fiber reflects a balance between the viscous and elastic forces upon it. The value of the diameter influence 
both of them: 

• The curvature Kf of a fiber of length Lf submitted to a force q per unit length scales as qL 2 / (E I): 
El is the bending modulus with I = 7rZ? 4 /64. The force necessary to achieve a given curvature scales 
therefore as D 4 . 

• In the present case, q is a viscous stress associated to the local relative velocity q also depends on D 
through the confinement ratio D/H : Increasing D/H results indeed in a stronger blockage of the flow 
normal to the fiber and, thus, enhances the normal viscous stress. This confinement effect is taken 
into account by introducing a factor A /f(D/H) in the proportionality relation between the viscous 
stress and both the relative velocity of the fiber and the fluid and the viscosity of the latter. The 
influence of D/H on the force on fixed circular cylinders in a similar flow geometry has been studied 
in Refs. 1261 - 1281 For the ratios D/H = 0.22 and 0.41 corresponding to D = 0.49 and 0.91mm, these 
studies provide respective values A^ = 30 and 73. 

As suggested in the previous section, the viscous forces are proportional to the local relative velocities 
of the different parts of the cylinders and the fluid in the regions where the streamlines are curved. These 
relative velocities may be estimated by multiplying dv/dn by the local distance between the ends of the 
fiber (of curvature Kf) and a streamline of curvature k s : this distance is then of the order of (k s — Kf)L 2 so 
that we take: q oc A^ p, (k s — Kf)L 2 (dv/dn). The factor A^ takes into account the confinement of the flow 
Using this expression in the estimation of the curvature Kf leads to: 


K f 


K s - Kf 


= cs p = c 


A /; [j, L* (dv/dn) 


El 


( 2 ) 


The right hand side of the equation is the Sperm number S p defined by Eq. ([!]) with the additional factor 
Ap • 

Fig. [8] shows the variation of Kf/n s as a function of S p . This time, the data corresponding to the two 
different fiber diameters collapse onto a same global trend. At low values of S p -C 100, the ratio Hf/n s 
increases steeply with S p (k//k s — S p /110 for the dashed line). At large values of S p > 500, Kf / k s becomes 
of the order of 1 and the profile of the fiber follows the streamlines. The dashed line corresponds to the best 
fit of these experimental data by Eq. 0 which is obtained for C = 1/110 so that Kf /k s = ^/(llO + S p ): 
the fitted curve follows well the trend of the experimental data. 

Two of the data points plotted in Fig. [5] correspond to the lower viscosity = 1 mPa.s. Like in Fig. 0 
(and for the same reasons), the corresponding values of Kf /k s are significantly lower. 


E. Influence of the length of the fiber 

All the experiments described above were performed with fibers of length Lf = 30 mm. Further mea¬ 
surements of the curvature and longitudinal velocity have been performed both for shorter and longer fibers 
(Fig. O: increasing the length Lf increases indeed the normal viscous stress (oc L 2 ) while it reduces the 
elastic force needed for a given bending (oc LJ 2 ). In agreement with Eq. 0, Kf/K s should then scale like 
L A f while it is -C 1. 

Fig. [9] displays the shapes of fibers of different lengths carried over the obstacle by a flow of given velocity 
and viscosity. While the curvature decreases clearly as Lf drops from 30 to 15 mm, its value for the longest 
fiber (Lf = 60 mm) is similar to that for Lf = 30 mm. As seen in Fig. [Gl bending is indeed localized to the 
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FIG. 8. Dimensionless curvature, Kf / k s of fibers of constant length: Lf = 30 mm as a function of the Sperm number 
S p = \^p,(dv/dn)L)/EI. Fiber diameter: D = 0.49 mm (□, ^), D = 0.91mm (o); fluid viscosity: p = 15mPa.s (o, 
□ ), n = 1 mPa.s (^). Dashed line: predictions from Eq. (j2j with C = 1/110. 



FIG. 9. Superimposition of fibers of different lengths Lf = 15 (dotted lines), 30 (continuous lines) and 60 mm 
(dashed lines) over two streamlines (thin grey lines) at two locations upstream (left) and downstream (right) of the 
step. Fluid velocity: U = 230mm.s _1 ; transverse locations of injection point: yi n i e t = 4.5 and 21.5 mm 


regions of maximum curvature of the streamlines. Therefore, for Lf = 60 mm, only a fraction of the length 
of the fiber is exposed to a large normal viscous stress: this exposed length is of the order of the width 
AW ~ 30 mm of each extremum of dv/dn in the curves of Fig. [3ju and is independent of Lf. For long fibers 
(Lf > AW), the relevant length to be used in the expression of S. p is therefore AW instead of Lf. For this 
reason, we replace the expression of S p from Eq. © by the more general expression: 

0 Ap /i£ 4 (dv/dn) /oN 

“ El ’ ( ) 

in which £ = Lf for Lf < AW and £ = AW for Lf > AW. The validity of these assumptions is tested in 
Fig. [TOl displaying the variation of Hf/n s with S p for fibers of different lengths and diameters. For the two 
shortest fibers Lf = 15 and Lf = 30 mm, one has Lf < AW so that both definitions of S p are valid; for 
Lf = 60 mm, a better collapse is clearly obtained by using the definition of Sp from Eq. © (A) than from 

Eq. ©. 

IV. DISCUSSION AND CONCLUSION 

The present experiments identified the parameters controlling the motion and deformation of a flexible 
fiber in a confined Hele-Shaw flow where a step-like obstacle induces variations of the magnitude and 
orientation of the velocity. 

Both the trajectory of the center of mass of the fiber and its orientation follow globally those of the 
corresponding streamline. This contrasts with pinned flexible fibers or fibers placed in a Poiseuille flow, 
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FIG. 10. Dimensionless curvature Kf/Ks as a function of S p for different lengths of the fibers: Lf = 15 mm ((0) 
D = 0.49mm), Lf = 30mm ((O) D = 0.49mm, (□) D = 0.91mm) and Lf = 60 mm ( D — 0.91mm) with AW 
(A) and Lf (a) as the characteristic length in the definition of Sp. Fluid viscosity: p = 15mPa.s. Dashed line: 
predictions from Eq. © with a proportionality factor 1/110. 


for which the fibers may cross the streamlines. Moreover, the velocity of the fibers is proportional to the 
fluid velocity at their center of mass with a constant only depending of the geometrical confinement. The 
variations of the curvature Kf of the fiber also follow closely those n s of the streamlines but, generally, with 
a smaller amplitude. 

Quantitatively, a key result of the present study is that the normalized curvature k//k s of the fibers 
is only a function of a combination S p (the Sperm number) of different parameters of the flow given by 
Eq. ©). Physically, relation © between Kf/n s and S p reflects the balance between viscous forces applied by 
the flowing fluid on the fiber and elastic forces opposing the resulting deformations. The viscous forces are 
proportional to the transverse velocity gradient, to the local distance between the fiber and the streamline 
and to the fluid viscosity. The elastic restoring forces, per unit length, are proportional to the curvature of 
the fiber, to its bending stiffness El oc D 4 and vary as £ -2 with its effective length £ for a given curvature: 
£ coincides with the length Lf only for short libers. For longer ones, £ must be taken equal in the expression 
of S p to a characteristic length of the variation of the flow. Here, if Lf > AW, £ is equal to the width AW 
of the region where the component dv/dn of the gradient of the velocity is large (Fig.©)); parts of the fibers 
outside this region do not influence bending. 

In contrast with many previous studies of the dynamics of flexible fibers, one deals here with a confined 
configuration of a fiber between two parallel walls: this keeps the motion of the fiber bi-dimensional and 
enhances the viscous forces due to the partial blockage of the flow: this effect is taken care of by introducing 
a dimensionless parameter depending on the diameter. The diameter influences however bending mainly 
through the D A factor in El which increases the rigidity considerably and reduces bending for fibers of 
larger diameter. 

The present study can be applied to the transport of flexible fibers in fractures containing obstacles 
occupying a part of the area: the quantitative information we obtained can be transposed to estimate the 
influence of the length and diameter of the fibers, of the fracture gap and the size in the fracture plane of 
the obstacles (modeled here as step reducing the width of an Hele-Shaw cell). 

In the present flow geometry, the velocity of the flow increases with distance which keeps the fiber under 
tension, except in localized regions. It will be important to investigate the opposite case in which buckling 
may, for instance, occur. The discussion presented above is largely based on a local balance between elastic 
and viscous forces. Non local effects might however occur for fibers submitted to a sequence of variations 
of curvatures of the streamlines over distances smaller than their length. Interactions between the different 
parts of the fiber may then be significant: this may have an important influence on transport in porous 
media or micro-fluidic circuits. 
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